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Signatures of tiie vulcanization transition — amorphous solidification induced by the random 
crosslinking of macromolecules — include the random localization of a fraction of the particles and 
the emergence of a nonzero static shear modulus. A semi-microscopic statistical-mechanical theory is 
presented of the latter signature that accounts for both thermal fluctuations and quenched disorder. 
It is found (i) that the shear modulus grows continuously from zero at the transition, and does so 
with the classical exponent, i.e., with the third power of the excess cross-link density and, quite 
surprisingly, (ii) that near the transition the external stresses do not spoil the spherical symmetry 
of the localization clouds of the particles. 
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Introduction: When a sufficient density of randomly lo- 
cated cross-links is imposed on a system of flexible linear 
macromolecules, an equilibrium phase transition (known 
as the vulcanization transition) occurs. At this transi- 
tion a liquid state is replaced by an amorphous solid 
state. This transition has two main equilibrium signa- 
tures: (i) a nonzero fraction of the monomers become 
localized around random mean positions and with ran- 
dom localization lengths (structure); and (ii) the system, 
as a whole, acquires a nonzero static shear modulus (re- 
sponse). The former signature has been discussed pre- 
viously; the purpose of the present Letter is to address 
the latter signature. Specifically, our aim is to develop a 
statistical-mechanical theory of the elastic properties of 
the amorphous solid state in the vicinity of the vulcaniza- 
tion transition. This theory incorporates both annealed 
(i.e. thermally equilibrating) and quenched random (i.e. 
cross-link specifying) variables. Its primary conclusions 
are: (a) that the amorphous solid [in the sense of signa- 
ture (i)] state emerging at the vulcanization transition 
is indeed a sohd [in the sense of signature (ii)]; (b) that 
the shear modulus vanishes continuously as the transition 
is approached, and does so with the third power of the 
excess cross-link density (i.e. the amount by which the 
cross-link density exceeds its critical value); and (c) that 
the shearing of the container associated with elastic de- 
formations does not lead to a shearing of the probability 
clouds associated with the thermal fluctuations of local- 
ized particles about their mean positions. 

The elastic properties of vulcanized matter and re- 
lated chemically-bonded systems, especially those near 
the amorphous solidification transition, have received 
considerable attention to date. Notable approaches in- 
clude the classical ones in which it was argued that 
near the transition the elastic entropy in the solid phase 
(and consequently the static shear modulus E) grow as 
the third power of the excess cross-link density e, i.e., 
i? ~ e* with i = 3. Subsequently, it was proposed that 
the amorphous solidification transition of polymer sys- 
tems be identified with a percolation process [^,||. Thus, 
the exponent t was identified with the critical exponent fi 



for percolation of conductivity (with ~ 2.0 in 3 spatial 
dimensions) . Subsequently it was observed that the elas- 
ticity percolation exponent for a random network is sub- 
stantially higher than /j, when the forces are central jj]. 

More microscopically oriented approaches to the elas- 
tic properties of vulcanized matter have also been made, 
in which macromolecular degrees of freedom feature ex- 
plicitly. Among these are the "phantom network" |^ 
and "affine network" [^ approaches, as well as the com- 
prehensive discussion of rubber elasticity by Deam and 
Edwards 0, and others ||]. These approaches focus on 
the well-cross-linked regime rather than the lightly-cross- 
linked regime near the vulcanization transition ^ . 

Experimentally, the exponent t has been addressed 
for several systems (although mostly for gelation rather 
than vulcanization): the results vary from t ^ 2 [^ to 
t > 3 [|ll|. This wide discrepancy is not understood. 

Stimulating though they certainly are, it must be rec- 
ognized that neither the classical [|l^,^,0 nor the perco- 
lation [^,^ approaches to the physics of vulcanized mat- 
ter explicitly include both crucial ingredients: thermal 
fluctuations and quenched disorder. In recent years, an 
approach to the vulcanization transition has been devel- 
oped [|l^ [l^ that takes into account both of these in- 
gredients in the context of a semi-microscopic model for 
flexible, randomly cross-linked macromolecules. This ap- 
proach is very much inspired by the work of Edwards and 
collaborators , as well as by concepts from the field 
of spin glasses. Emerging from this more recent approach 
has been a detailed picture of the structure of the amor- 
phous solid state near to the vulcanization transition, 
including, in particular, an explicit form for the distri- 
bution of localization lengths. What has not yet been 
elucidated using this approach is the second signature 
of the vulcanization transition, namely the emergence of 
static response to shear deformations. This issue is the 
focus of the present Letter. 

Model: At the heart of the theory of the structure of 
the amorphous solid state [|l^ is the analysis, employ- 
ing the techniques of replica statistical mechanics, of a 
semi-microscopic model of A'' macromolecules subject to 
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random cross-linking constraints. This analysis leads to 
an order parameter appropriate for diagnosing the 
amorphous solid state, as well as a Landau free energy 
J-n{{^^.}) in terms of this order parameter. A detailed 
review of this theory has been given in Ref. The 
order parameter is defined via 



(1) 



Here, hatted vectors denote replicated collections of vec- 
tors, viz., V = {v°, v^, • • ■ , v"}, their scalar product be- 
ing V ■ w = J2a=o ■ ^^'^ trajectories {c(s)}^j 
are the semi-microscopic configurations of the replicated 
macromolecules (where < s < 1 is the arclength in 
units of the total arclength). {■ ■ denotes an av- 

erage for an effective pure (i.e. disorder-free) system of 
n + 1 coupled replicas of the original system. To model 
the disorder we make the Deam-Edwards assumption |^] 
that the statistics of the cross-links is determined by the 
instantaneous correlations of the uncrosslinked system. 
This leads to the need to work with the n — > limit of 
systems of n -|- 1 (as opposed to n) rephcas. The addi- 
tional replica, labeled by a = 0, represents the degrees 
of freedom of the original system before cross-linking, or, 
equivalently, describes the cross-link distribution. Con- 
sequently, any external strain applied to the system after 
the permanent constraints have been created will affect 
replicas a — 1, . . . , n, but not replica a = ||^. Thus, 
the order parameter measures the correlations between 
the positions of individual particles before and after the 
deformation is applied. 

In the saddle-point approximation ||l9|, the disorder- 
averaged free energy / (per particle) in a d-diniensional 
system is obtained by minimizing the replicated free- 
energy functional ^„({^2j,}) ||20|| : 



/ = d lim min :r„({17-}). 



(2) 



As discussed in detail in Ref. |]T^, the minimization in 
Eq. (H) yields the liquid-amorphous-solid phase transi- 
tion at a certain critical value of the cross-link density. 
We parametrize the excess cross-link density beyond this 
critical value by the control parameter e. As the tran- 
sition is continuous (i.e. near the critical point the gel 
fraction is small and the typical localization length of lo- 
calized particles is large), ^„({r2^}) can be expanded in 
powers of the order parameter and gradients, with only 
low orders needing to be retained pi[ : 

The symbol ^ denotes a sum over replicated wave- 
vectors that contain at least two nonzero component- 
vectors k" The saddle-point equation for the free- 
energy functional near the transition is exactly solved by 
the following hypothesis : 



W{k) EE / dTp{T)e-^'/^^ . (4b) 
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The physical motivation for this hypothesis comes from a 
picture in which a fraction q of the monomers are local- 
ized around random mean positions hi{s) about which 
they execute harmonic thermal fluctuations over random 
localization lengths ^i(s). Furthermore, the mean posi- 
tions are assumed to he homogeneously distributed over 
the sample, and the localization lengths are character- 
ized by the statistical distribution 2^~^p(^~^). Thus, 
delocalized and localized particles are, respectively, rep- 
resented by the first and second terms on the RHS of 
Eq. (^). The (5- factor in the second term comes from 
the homogeneity of the distribution of mean positions. 
The function W^{k), which we refer to as the continuous 
part of the order parameter, encodes all the information 
about thermal fluctuations (the superscript u standing 
for "unstrained"). The hypothesis ( [4a|) and (4b) satisfies 
the saddle-point equations provided that ^6 17 



0=-2qe + 3q\ 



(5a) 



= (1 - 61) 7r(6l) - / d9'7r{6')TT{9-9'), (5b) 
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where t:{9) is an e-independent scaling function such that 
p(^~^) = (2/e) 7r(2/e^^), and satisfies the boundary con- 
dition j^d9'K{9) = 1. Equation ( ^a[ ) determines the lo- 
calized fraction q: for e < we obtain q — (i.e. the 
liquid phase, which has a vanishing static shear modu- 
lus); for e > we obtain q = 2e/3, corresponding to the 
amorphous solid state, which is the state on which we 
shall focus from now on. 

Response to shear strain: We now set about determining 
the free-energy cost associated with making static shear 
deformations of the system. To do this, we consider the 
effect of changing the shape of the container (on which we 
have imposed periodic boundary conditions). We char- 
acterize the deformation by the {d x d) matrix S, which 
describes the change in position of any point b at the 
boundary of the system as follows: b ^ S • b. For ex- 
ample, for d = 3 and for a deformation in which the x, 
y and z cartesian components of the position vector are, 
respectively, elongated by the factors A^^, Aj, and A^, the 
matrix S has the form diag(A2;, Aj,, Az). As we are con- 
cerned with the free-energy cost of pure shear strains, we 
shall assume that the deformation leaves the volume V 
of the system unchanged, i.e., Det S = 1. For considering 
infinitesimal strains, it is convenient to define the (sym- 
metric) strain tensor J = i(S -|- S""") — I. Here S"^ is the 
transpose of S, and I is the identity matrix. As Det S = 1 
we have Tr J = 0, to first order in the deformation. 

Before taking the thermodynamic limit, the system is 
finite in extent, so that the Fourier representation of any 
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function of position consists of a superposition of plane 
waves with discrete wave-vectors. In particular, the or- 
der parameter (which is a function on replicated Fourier 
space) is only defined at a discrete set of points. Now, 
under strain the boundaries in position space are dis- 
placed and, as a consequence, the discretization in repli- 
cated Fourier space changes. As mentioned above, any 
external strain applied to the system after the perma- 
nent constraints have been created will affect replicas 
a = 1, . . . ,n, but not replica a = 0. Therefore, the 
change in the discretization of the wave vectors occurs 
only for a ~ 1, . . . ,7i, but not a — 0. For convenience, 
we shall use the symbols i?" and to denote the sets 
of allowed replicated wave- vectors corresponding, respec- 
tively, to the unstrained and strained systems. 

Conceptually, there are two sources for the change in 
free energy [Eq. (^] under deformation: the change in the 
expression for the free energy functional itself, and the 
consequent change in the value of the order parameter 
that solves the saddle-point equation. The free-energy 
functional for the strained system ^^({fi^}) is obtained 
by repeating, step-by-step, the construction of the free- 
energy functional for the unstrained system J>i({J7j,}). 
The result is that the coefficient in front of each 
term is unaltered, the only change being the replacement 
of each sum J2k °f discrete replicated 

wave-vectors (i.e. k € i?") by a sum over the new set of 
discrete replicated wavevectors J2k£R^- -^^ ^ result, the 
saddle-point equation for the strained system becomes 
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We now obtain the order parameter for the strained 
system by making a physically motivated hypothesis sim- 
ilar to the one made for the unstrained system. First, for 
each localized monomer in the unstrained system we en- 
vision that its old mean position hi{s) is displaced to a 
new mean position S • hi{s) + ri(s), where S • hi{s) is 
the afSne displacement of the old position [|| and ri(s) is 
a random additional displacement, which we take to be 
uncorrelated with hi{s). With the assumption that (as 
in the unstrained system) there is no correlation between 
the extent (including shape) of the thermal fluctuations 
of a monomer about its mean position and the mean po- 
sition itself, we arrive at the hypothesis 



where W^{k) is the continuous part of the order pa- 
rameter in the strained system. Now, to construct a 
form for W^{k) we consider a conjecture for the form 
of (e*'' '^'^*-')^ (i.e. the thermal expectation values of the 
Fourier-transformed individual particle densities in the 
strained (s) system for a specific disorder realization x)- 

exp(ik-{S-b,(s) -f r,(s)}) 

xexp(-e'(s)k-{I + 77.(s)J}.k/2). (8) 



We expect the gaussian probability cloud to be isotropic, 
except for a correction due to the distortion. For in- 
finitesimal distortions, this correction should be propor- 
tional to J and have a random magnitude ??i(s). For ex- 
ample, if rii(s) — 2 then the probability cloud would have 
been affinely distorted. On the other hand, if rii{s) — 
then the probability cloud would remain spherical. As- 
suming, further, that the random displacement ri{s) also 
has a probability distribution shaped by a combination 
of I and J, and expanding to lowest nontrivial order in 
the deformation p3, we obtain the hypothesis: 



W\k) = q dee 
Jo 



m{0) 



^k"-J-k"). (9) 



Here, w(9) is a second scaling function, which describes 
the change in the continuous part of the order parameter 
due to the deformation. 

Alternative motivation for the form of W^{k) runs as 
follows. Let us assume that for small strains W^{k) is 
unchanged under simultaneous rotations of S and k. As 
it is only a function of fc^ this property certainly holds 
for W^(k), and it therefore also holds for the difference 
between W^{k) and W^{k). To first order in J this dif- 
ference can only contain the following terms: (i) a linear 
function of Eq=i k"- J-k" and (ii) a product of an invari- 
ant linear function of J with an invariant function of k. 
The only quantity linear in J and invariant under rota- 
tions is Tr J, which is zero for infinitesimal shear strains, 
as mentioned above. Thus we recover Eq. (^. 

By inserting the hypothesis given by Eqs. (|^) and (0) 
into the saddle-point condition (^, we recover Eqs. ( |5a| ) 
and (Bq) for q and 7r(6'), together with the condition: 



6*2 dw 
Y~d9 



de' 0' m{0')Tr{e 



(10) 



The boundary condition limg-^oo 0'^zu{0) — stems from 
the fact that, by Eq. (|l|), hm|j,|^^ fl^, ~ 0. The only so- 
lution of Eq. (|o|) that satisfies the boundary condition is 
the null function: zu{0) = |2^. This result implies the 
first (and, a priori, the most surprising) result of this Let- 
ter: the continuous part of the order parameter does not 
change to first order in the strain, i.e., W^(k) = W^{k). 
This conclusion is consistent with the phantom network 
picture §J|. It also suggests that W^{k) = Wik) for 
finite (and not merely infinitesimal) deformations, and 
indeed the resulting order-parameter hypothesis turns 
out to satisfy the saddle-point equation for arbitrarily 
strained systems. 

We now have all the ingredients necessary to compute 
the change in the free energy, to leading order in e, due 
to the deformation of the system: 



-d lim 



2^ 
27 



-TrfS-S' 
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where n| and fl^ are, respectively, the saddle-point val- 
ues of the order parameter for the strained and un- 
strained systems. Thus we can extract the value of 
the static shear modulus for the amorphous solid state 
near the solidification transition (with physical units re- 
stored): E = k-QTNCe", where k-B is Boltzmann's con- 
stant, T is the temperature, and C is a model-dependent 
positive constant. Hence we see that the static shear 
modulus near the vulcanization transition is character- 
ized by the exponent t = 3, in agreement with the classi- 
cal result A simple scaling argument, viz., that the 
modulus should scale as two powers of the order param- 
eter (q^) and two powers of the gradient (Ctyp); leads to 
the same value for t. 

Concluding remarks: We have presented a microscopic 
derivation of the static elastic response of a system of ran- 
domly cross-linked macromolecules near the amorphous 
solidification transition. In the picture that emerges, it 
is seen: (i) that the amorphous solid state, which was 
previously shown to be characterized structurally by the 
localization of a nonzero fraction of particles, is also 
characterized by having a nonzero static shear modu- 
lus; (ii) that the static shear modulus scales as the third 
power of the excess cross-link density (beyond its value 
at the transition) |2^] ; and (iii) that the form of localiza- 
tion exhibited by the particles is left unchanged by the 
strain. It is, however, not implausible that strain-induced 
changes would emerge from a more detailed analysis of 
the effects of the excluded- volume interaction, at least at 
higher crosslink densities. Being dependent only on the 
form of the free-energy functional [^,^, and not any 
specific semi-microscopic model, the approach to elas- 
ticity described here should be generally applicable not 
only to systems of randomly cross-linked flexible macro- 
molecules, but also to other equilibrium amorphous solid 
forming systems. 
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